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Gravitinos are very promising candidates for the cold dark matter of the Universe. Interestingly, to 
achieve a sufficiently long gravitino lifetime, -R-parity conservation is not required, thus preventing 
any dangerous cosmological influence of the next-to-lightest supersymmetric particle. When R- 
parity is violated, gravitinos decay into photons and other particles with a lifetime much longer 
than the age of the Universe, producing a diffuse gamma ray flux with a characteristic spectrum 
that could be measured in future experiments, like GLAST or AMS-02. In this letter we compute 
the energy spectrum of photons from gravitino decay and discuss its main qualitative features. 
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There is mounting evidence that dark matter is ubiq- 
uitous in our Universe [ij. Since the necessity of dark 
matter wa.s realized, many particle physics candidates 
have been proposed. Among the most interesting stands 
the gravitino Q, the supersymmetric counterpart of the 
graviton, which arises when global supersymmetry is pro- 
moted to a local symmetry. If the gravitino is the lightest 
supersymmetric particle (LSP), it constitutes an excel- 
lent candidate for the cold dark matter of the Universe. 
The thermal relic density of the gravitino is calculable in 
terms of very few parameters, the result being [sl 




(1) 

while the relic density inferred by WMAP for the ACDM 
model is ^cBuh^ ~ 0.1 In this formula, Tr is the 
reheating temperature of the Universe, 7713/2 is the grav- 
itino mass and rrig is the gluino mass. It is indeed re- 
markable that the correct relic density can be obtained 
for typical supersymmetric parameters, 7713/2 ~ 100 GeV, 
77ig 1 TeV, and a high reheating temperature, Tr ~ 
10^° GeV, as required by baryogenesis through the mech- 
anism of thermal leptogenesis 

Whereas the gravitino as the LSP leads to a cos- 
mology consistent with observations, the cosmology of 
the Next-to-Lightest Supersymmetric Particle (NLSP) is 
much more problematic. In supersymmetric model build- 
ing, in order to prevent too rapid proton decay, it is 
common to invoke a discrete symmetry called i?-parity. 
When i?-parity is exactly conserved, the NLSP can only 
decay into gravitinos and Standard Model particles with 
a decay rate strongly suppressed by the Planck mass. As 
a result, the NLSP is typically present in the Universe 
at the time of Big Bang nucleosynthesis, jeopardizing 
the successful predictions of the standard nucleosynthesis 
scenario. In most supersymmetric scenarios, the NLSP is 
either a neutralino or a stau. On one hand, if the NLSP 
is a neutralino, its late decay into hadrons can dissociate 
the primordial elements Q. On the other hand, if the 
NLSP is a stau, it can form a bound state with '*He, cat- 
alyzing the production of ^Li 0]. As a result, the abun- 
dance of ^Li is increased by a factor 300-600, in stark 
conflict with observations Isfl. 



Several scenarios have been proposed that circumvent 
the above-mentioned difficulties f9|. The simplest, albeit 
the most radical one, is based on the assumption that 
i?-parity is not exactly conserved In fact, although 
experiments set very stringent bounds on i?-parity vio- 
lation, there is no deep theoretical reason why it should 
be exactly conserved. If i?-parity is mildly violated, the 
NLSP decays into Standard Model particles well before 
the first nucleosynthesis reactions take place, thus not 
posing a jeopardy for the Standard Model predictions. 
Remarkably, even though the gravitino is no longer sta- 
ble when i?-parity is violated, it still constitutes a viable 
dark matter candidate To be precise, a consistent 

thermal history of the Universe with gravitino dark mat- 
ter, thermal leptogenesis and successful primordial nucle- 
osynthesis requires the lepton number violating Yukawa 
couplings to lie between 10"^'' and 10"'' which trans- 
lates into gravitino lifetimes in the range 10^'^ — 10^^ s 
for 7rt3/2 ~ 100 GeV, which are much longer than the age 
of the Universe. Nevertheless, gravitino decays could be 
happening at a sufficiently high rate for the decay prod- 
ucts to be detectable in future experiments. 



In this letter we will concentrate on the photons pro- 
duced in gravitino decays, although in general other sta- 
ble particles are also produced, such as electrons, protons, 
neutrinos and their antiparticles. Demanding a high re- 
heating temperature for the Universe as suggested by 
thermal leptogenesis, Tr > 10^ GeV [ll, it follows from 
Eq. Ill that the gravitino mass has to be 7^13/2 > 5 GeV 
for typical gluino massses. Consequently, we expect the 
photons from gravitino dark matter decay in the energy 
range of a few GeV, i.e. in the gamma ray energy range. 

The Energetic Gamma Ray Experiment Telescope 
(EGRET) aboard the Compton Gamma Ray Obser- 
vatory measured gamma rays in the energy range be- 
tween 30 MeV to 100 GeV. After subtracting the galac- 
tic foreground emission, the residual flux was found to 
be roughly isotropic and thus attributed to extragalac- 
tic sources. The first analysis of the EGRET data by 
Sreekumar et al. [l3| gave an extragalactic flux with an 
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energy spectrum described by the power law 
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in the energy range 50 MeV-10 GeV. The improved anal- 
ysis of the galactic foreground by Strong et al. opti- 
mized in order to reproduce the galactic emission, shows 
a power law behavior between 50 MeV-2 GeV, but a clear 
excess between 2-10 GeV, roughly the same energy range 
where one would expect a signal from gravitino decay. 
Although it is very tempting to look for explanations for 
this excess in terms of gravitino decays, in view of all the 
systematic uncertainties involved in the extraction of the 
signal from the galactic foreground, we will not attempt 
to fit our predicted flux to the EGRET data. Nonethe- 
less, we will show later the EGRET data superimposed 
with our predicted flux for comparison. 

The total gamma ray flux received from gravitino dark 
matter decay receives two main contributions. The first 
one stems from the decay of gravitinos at cosmological 
distances, giving rise to a perfectly isotropic extragalactic 
diffuse gamma ray background. Defining dN^/dE as the 
gamma ray spectrum produced in the gravitino decay, 
the flux received at the Earth with extragalactic origin 
has the following expression: 
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z, z being the redshift, and 
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Here, ^^3/2, and flA are the gravitino, matter and 
cosmological constant density parameters, respectively, 
Pc is the critical density, r3/2 the gravitino lifetime, and 
Hq the present value of the Hubble parameter. 

In addition to the cosmological contribution, the to- 
tal gamma ray flux also receives a contribution from the 
decay of gravitinos in the Milky Way halo. This contri- 
bution reads: 
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The integration extends over the line of sight, so the 
halo contribution has an angular dependence on the di- 
rection of observation, yielding an anisotropic gamma 
ray flux that in EGRET could resemble an isotropic 
extragalactic flux. Namely, in the energy range 0.1-10 
GeV, the anisotropy between the Inner Galaxy region 
(|6| > 10°, 270° < I < 90°) and the Outer Galaxy region 
(\b\ > 10°, 90° < I < 270^) is just a 6%, well compatible 
with the EGRET data 1J| (for a detailed discussion ap- 
plied to GLAST, see In Eq. |[5|), Phaio stands for the 
dark matter distribution in the Milky Way halo. For our 



numerical analysis we will adopt a Navarro- Frenk- White 
density profile [1^ 
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where r is the distance to the Galactic center, Tc — 20 kpc 
is the critical radius and ph — 0.33 GeVcm"'^. 

In Eqs. I|3|5p the only undetermined quantity is the 
energy spectrum of photons produced in the gravitino 
decay, dN^/dE, which depends crucially on the gravitino 
mass. If the gravitino is lighter than the bosons, it 
decays mainly into a photon and a neutrino by means of 
the photino-neutrino mixing that arises when i?-parity is 
violated [lH. Therefore, the spectrum is simply 



dN^ 
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For this case, it was found in [lOj, ll5| that the total 
gamma ray fiux received is dominated by the monochro- 
matic line from the decay of gravitinos in our Milky Way 
halo, while the redshifted line from the decay of graviti- 
nos at cosmological distances is somewhat fainter. 

On the other hand, if the gravitino is heavier than the 
or bosons, new decay modes are open. In ad- 
dition to the decay mode into a photon and a neutrino 
that follows from the photino-neutrino mixing, U^,^, the 
gravitino can also decay into a boson and a charged 
lepton, through the mixing charged wino-charged lepton, 
U^,g, or into a Z° boson and a neutrino, through the mix- 
ing zino-neutrino, C/^^. The decay rates can be straight- 
forwardly computed from the interaction Lagrangian of 
a gravitino with a gauge boson and a fermion [lj|. Ne- 
glecting the masses of the final fermions, the result for 
each decay mode can be approximated by: 
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where f{x) — 1 , g„ 

The fragmentation of the and the Z° gauge bosons 
will eventually produce photons, mainly from the decay 
of neutral pions. We have simulated the fragmentation 
of t he g auge bosons with the event generator PYTHIA 
6.4 and calculated the spectra of photons in the 
and Z° channels, which we denote by dN^ /dE and 
dN^ /dE, respectively. The total spectrum is given by: 
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The branching ratios in the diflerent decay channels 
are determined by the size of the i?-parity breaking mix- 
ing parameters, U^^, U^^ and Ut^^, and by the kine- 
matical function f{x) defined after Eq. ([8]). The mixing 
parameters stem from the diagonaHzation of the 7x7 
neutraUno-neutrino and 5x5 chargino-charged lepton 
mass matrices, whose expHcit form can be found in the 
vast existing hterature on i?-parity violation [l^. The 
precise expression for the mixing parameters in terms 
of the i?-parity breaking coupHngs in the Lagrangian is 
fairly cumbersome and will not be reproduced here. How- 
ever, to derive the branching ratios, only the ratio among 
them is relevant, and not their overall value. 

To derive the relation between U^^ and U^^, we first 
note that the photino does not couple directly to the neu- 
trino (since neutrinos do not couple to photons). Never- 
theless, an effective photino-neutrino mixing is generated 
through the mixing photino-zino and the mixing zino- 
neutrino. The result reads: 
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Therefore, the relation between U^i, and U follows from 
the 2x2 gaugino sub-block of the neutralino mass matrix, 
that in the (— 17, —iZ) basis reads 
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Here, Mi and M2 are the U{1)y and SU{2)l gaugino 
masses, and Cw {sw) denotes the cosine (sine) of the 
weak mixing angle. Therefore, 
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which depends only on the gaugino masses at low ener- 
gies. Assuming gaugino mass universality at the Grand 
Unified Scale, Mx = 2 x lO^^GeV, we obtain at the 



electroweak scale M2/M1 

0.31IC/; 



1.9, which yields \U^ 



The mixing parameter C/f^^, on the other hand, is re- 
lated to J7^^ by SU{2)l gauge invariance. The relation 
approximately reads: 
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M2 is the wino mass at the electroweak 



scale. Using Eq. ifTTI) . we finally obtain 
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which under the assumption of gaugino mass universal- 
ity at Mx yields \Uyy^\ ~ 1.09|J7^i/l- Hence, under this 
assumption, the mixing parameters are in the ratio 



wei 



1 : 3.2 : 3.5 



TABLE I: Branching ratios for gravitino decay in different 
-R-parity violating channels for different gravitino masses. 

7713/2 BR(V;3/2 ^v) BR(V;3/2 ^ Wt) BR(V;3/2 ^ Z%) 



10 GeV 
85 GeV 
100 GeV 
150 GeV 
250 GeV 



1 

0.66 
0.16 
0.05 
0.03 





0.34 
0.76 
0.71 
0.69 






0.08 
0.24 
0.28 



and thus the branching ratios for the different decay 
modes only depend on the gravitino mass (see Table H). 

Once the spectrum of photons from gravitino decay has 
been computed, Eq. |[9]), it is straightforward to calculate 
the gamma-ray fiux received at the Earth from our local 
halo and from cosmological distances, by using Eqs. (|3|5p . 
Assuming universality of gaugino masses at high energies, 
the photon fiux received from gravitino decay depends 
essentially on the gravitino mass, which determines the 
shape of the energy spectrum, and the gravitino lifetime, 
which determines its overall normalization. 

In Fig. [T] we show the different contributions to the 
gamma ray fiux for 7713/2 = 150 GeV and T3/2 ~ 2 x lO^^s. 
To compare our results with the EGRET data [l3|, also 
shown in the figure, we have averaged the halo signal 
over the whole sky excluding a band of ±10° around 
the Galactic disk, and we have used an energy resolu- 
tion of 15%, as quoted by the EGRET collaboration 
in this energy range. The energy resolution of the de- 
tector is particularly important to determine the width 
and the height of the monochromatic line stemming from 
the two body decay -03/2 1^. The three contribu- 
tions are dominated by the halo component, the ex- 
tragalactic component being smaller by a factor of 2- 
3. Finally, we have adopted an energy spectrum for 
the extragalactic background described by the power law 
= 4 X 10-^ (Jv)"°' (cm^^tr s)-iGeV, in or- 
der to provide a qualitatively good agreement of the total 
fiux received with the data. 

The predicted energy spectrum shows two qualitatively 
different features. At energies between 1-10 GeV, we 
expect a continuous spectrum of photons coming from 
the fragmentation of the gauge bosons. As a result, the 

(13) predicted spectrum shows a departure from the power 
law in this energy range that might be part of the ap- 
parent excess inferred from the EGRET data by Strong 
et al. [l3|. The upcoming space-based gamma ray ex- 
periments GLAST and AMS-02 will measure the energy 
spectrum with unprecedented accuracy, providing very 

(14) valuable information for the scenario of decaying grav- 
itino dark matter. 

In addition to the continuous component, the energy 
spectrum shows a relatively intense monochromatic line 
at higher energies arising from the decay channel tp3/2 
^v. This line could be observed not only by GLAST 

(15) or AMS-02 in the diffuse gamma background, but also 
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FIG. 1: Contributions to the total gamma ray flux for m3/2 = 
150 GeV and rg/a ~ 2 x 10^*^8 compared to the EGRET data. 
In dotted lines we show the photon flux from the fragmenta- 
tion of the Z boson, in dashed lines from the fragmentation 
of the W boson, and in dot-dashed lines from the two body 
decay tp3/2 — > 7!^. The background is shown as a long dashed 
line, while the total flux received is shown as a thick solid line. 



by ground-based Cherenkov telescopes such as MAGIC 
(with an energy threshold of 70 GeV) or VERITAS (50 
GeV) in galaxies such as M31 [l5| . 

The intense gamma Hne is very characteristic of this 
scenario, and the observation of this feature with the 
right intensity would support the gravitino dark mat- 
ter decay hypothesis. While scenarios with neutralino 



dark matter also predict a continuous spectrum and a 
monochromatic line coming from the annihilation chan- 
nels x^X^ ~' 111 ^1 01, these channels only arise at one 
loop level, and thus the intensity of the monochromatic 
line is greatly suppressed compared to the continuum. 
One should note, however, that the presence of an intense 
gamma Hne is not unique to the scenario with decaying 
gravitino dark matter and is also expected, for example, 
from the annihilation of inert Higgs dark matter [2l[. 

To summarize, in this letter we have computed the 
gamma ray flux from gravitino dark matter decay in sce- 
narios with i?-parity violation. These scenarios are very 
appeahng theoretically, as they naturally lead to a his- 
tory of the Universe consistent with thermal leptogenesis 
and primordial nucleosynthesis. The predicted flux es- 
sentially depends on two parameters: the gravitino mass, 
which determines the shape of the energy spectrum, and 
the gravitino lifetime, which determines its overall nor- 
maHzation. If the gravitino is Hghter than the and 
gauge bosons, the predicted energy spectrum is essen- 
tially monochromatic. On the other hand, if it is heavier, 
the energy spectrum consists of a continuous component 
and a relatively intense gamma ray line. This gamma ray 
flux might have already been observed by EGRET. Fu- 
ture experiments, such as GLAST, AMS-02 or Cherenkov 
telescopes, will provide unique opportunities to test the 
decaying gravitino dark matter scenario. 
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